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Abstract 

Modified dispersion relations from effective field theory are shown to alter the Chandrasekhar 
mass limit. At exceptionally high densities, the modifications affect the pressure of a degenerate 
electron gas and can increase or decrease the mass limit, depending on the sign of the modifica- 
tions. These changes to the mass limit are unlikely to be relevant for the astrophysics of white 
dwarf or neutron stars due to well-known dynamical instabilities that occur at lower densities. 
Generalizations to frameworks other than effective field theory are discussed. 
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I. INTRODUCTION 



The principle of Lorentz Invariance is at the heart of the contemporary formulation of 
physical theory, in particular the Standard Model and general relativity. Given the central- 
ity of Lorentz Invariance (LI) it is wise to explore as many avenues as possible that test 
this principle. One avenue investigated in recent years is the physics of modified dispersion 
relations (MDR). A remarkable result of these studies is that astrophysical data significantly 
limit Planck scale - suppressed modifications (see e.g. for reviews). The successful 

limits on the modifications are due to the high degree of sensitivity of particle process thresh- 
olds to Lorentz Violation (LV). Given the delicate interplay between the modifications of the 
particles involved, it is helpful to look for systems in which the effects of the modification 
on a single particle type are isolated. An apparently ideal system is the physics of a degen- 
erate electron gas since the pressure of such a gas supports the gravitational attraction of 
white dwarfs. Further, the Chandrasekhar mass limit, about 1.4M , of white dwarf stars is 
obtained in the ultra- relativistic limit, precisely where modifications of dispersion relations 
are expected to be large. 

Modified dispersion relations often take the form of an expansion in LV terms 

where the parameters Ki may differ for different particle species and Mp is the Planck mass, 



which we take to be M P = ^h/(A7iG) « 3.45 x 10 27 eV (c = 1). In such MDR models 
the usual energy-momentum conservation laws hold so there is a preferred frame, which we 
take to be the one where the cosmic microwave background is isotropic 181]. While these 
modified dispersion relations could be viewed simply as a phenomenological expansion to 
test LV, this form has been suggested in a variety of settings including string theory tensor 
vacuum expectation values, heuristic calculations of the semiclassical limit of loop quantum 
gravity, spacetime foam, non-commutative geometry, analogs of emergent gravity, and some 
braneworld models [l, 2|. 

Given the energy scale of the modifications, it might seem that testing such modifications 
might simply be impossible, However even 

process thresholds are highly sensitive to these modifications. Using several particle pro- 
cesses and observed energies, much of the parameter space is ruled out 0, Q] . For instance 
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the recent work of Macione et. al. achieves a limit on the parameter space of electrons of 
less than 10" 5 Q. 

In this paper we report on numerical solutions to the exact equations for the Chan 



drasekhar mass limit with modified dispersion relations. We extend the analysis of 
relaxing the unphysical assumption of constant density. We find significant differences with 
the reported results: The mass limit may be raised or lowered depending on the sign of 
the modifications in the electron dispersion relation; and physical equilibrium radii exist for 
both signs of the MDR parameter k. Finally we show that, despite the existence of new mass 
limits, there would be no effect on white dwarf astrophysics. The effects are only important 
in the Planck-scale regime far above astrophysically accessible densities. 

This paper is organized as follows. The next two sub-sections are devoted to discussions 
of MDR and the mass limit. In section [XT] we derive the corrections to the mass limit due 
to MDR. In the final section II I II we summarize the results, compare with results previously 
reported in 6|, and comment on the applicability of the calculation to other frameworks 
with modified dispersion relations. 

A. Modified Dispersion Relations 

To achieve precise limits on the parameters it is necessary to have some additional knowl- 
edge of the dynamics of the field theories. This may be achieved in the context of effective 
field theory, where effective field theory is used to determine the (non-renormalizable) mass 
dimension five (or higher) LV operators. 

Myers and Pospelov found that there are essentially only three operators that simulta- 
neously break local LI and preserve gauge and rotation invariance. Introducing a preferred 
frame field n a , the mass dimension 5 operators are j3] 

C -n a F ac n ■ d(n b F bc ) + _L n > 7 a(6 + M{n • <9)V (2) 



M P v " ' 2M P 

where (, & are dimensionless parameters, F ab is the dual of the usual F ab = d a A b — d b A a of 
Maxwell theory. The dispersion relations for fermions become 



E±=p 2 + m 2 + r ]± ^- (3) 

IVl p 

where the different signs correspond to the two helicity states. Typically, the modifications 
become important when they become comparable to the mass terms. For an electron this 



occurs at an energy scale of yJm^Mp ~ 10 TeV - significantly below the Planck scale. 
While these energy scales are beyond terrestrial experiments, astrophysical processes at 
these energy scales effectively constrain the cubic modifications 3|-[5j]. 

The fact that electron cubic MDR have two parameters leads to new effects. At sufficiently 
high energies electrons are unstable to helicity decay jl]. If r\- > r] + then the negative helicity 
electrons will decay producing a photon and positive helicity electron. While there is no 
kinematic threshold for this process, an effective threshold may be derived by studying 
the reaction rate. The resulting effective threshold is (m iMp/Ar]) 1 / 3 or about 10 TeV for 
Arj = \rj + — r]_\ < 0.4, as determined by photon stability [1]. The lifetimes of the negative 
helicity electron decrease around this effective energy; lifetimes for a 1 TeV electron is about 
1 s while for a 50 TeV electron it is reduced to about 10~ 9 s However the lifetimes increase 
at higher energies due to the fact that the states become more chiral. The MDR effects on 
the Chandrasekhar mass only occur at such high energies that the electron population is 
effectively in one helicity state. Assuming that the astrophysical processes, e.g. accretion, 
driving the increase in electron energy are sufficiently slow as compared to the lifetimes of 
the helicity states all the electrons would have positive helicity. Since the same analysis 
applies for the other case, rj + > r/_, the resulting degenerate electron gas will all be in the 
helicity state associated with the smallest parameter. Hence we will denote the electron 
parameter with k, where k = mm(r] + ,i]^). 

For the purposes of studying the effect on the physics of compact stars we assume k is 
order one and write the MDR as 



E 2 = p 



1 + K 



V 

Mr 



+ m. (4) 



While the case of dimension 5 operators and cubic modifications (n = 1) is the most interest- 
ing, since the analysis easily generalizes we present the mass limit calculation for arbitrary 
n. The range of validity for this effective description is p < Mp. 



B. Chandrasekhar Mass Limit 



White dwarf stars are an end stage of stellar evolution for stellar masses less than ~ 8M Q . 
As the star cools and contracts, the electron wavefunctions begin to overlap appreciably. As 
the star contracts further, the electrons are forced into higher and higher momentum states 
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to satisfy the Pauli exclusion principle, forming a degenerate electron gas with relativistic 
Fermi energy. Because the thermal energy is so much lower than the Fermi energy the 
temperature is effectively zero. The Chandrasekhar mass limit, about 1.4M , of white 
dwarfs stars arises when the Newtonian gravitational attraction due to the nuclei (often 
comprised of carbon and helium) is balanced by the outwards pressure of degenerate, cold 
electrons. The mass limit is obtained for an ultra-relativistic Fermi gas, when E/m e ~ 10 4 
with central densities in excess of p ~ 10 6 g cm" 3 . Thus the Chandrasekhar mass limit is 
derived in the framework of Newtonian gravitation and the statistical mechanics of an ideal, 
non-interacting, ultra-relativistic gas of electrons at T = 0. 

The Chandrasekhar mass limit is largely of theoretical interest. Other physics comes into 
play at high densities introducing processes known as dynamical instabilities. Depending 
on the details of the star's composition, two processes transform the physics of the star. 
At densities p ~ 10 9 — 10 11 the electrons acquire sufficient energy to induce inverse (5- 
decay and the star turns into a neutron star in a process called neutronization. At the 
critical density of ~ 10 10 g cm" 3 the onset of general relativistic instabilities causes the 
star to collapse. These threshold densities are composition dependent so, for instance, the 
general relativistic instability is irrelevant to iron white dwarfs, which undergo neutronization 
before gravitational collapse. For helium and carbon white dwarfs, the situation is reversed, 
instabilities due to general relativity occur before neutronization. 

II. CHANDRASEKHAR MASS LIMITS WITH MODIFIED DISPERSION RELA- 
TIONS 

To determine the actual mass limit - with or without modifications - a numerical analysis 
must be done. In calculating the electron degeneracy pressure, which supports the gravita- 
tional attraction, we assume that the star is isotropic and assume that the velocity of the 
star is not large in the preferred frame. The pressure is given by 



where n e is the electron number density. We use the group velocity in the MDR case, so 
that 




(5) 




(6) 
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As the temperature is effectively zero, all the states below the Fermi momentum, pp, are 
occupied. Hence the pressure in the continuum limit simply becomes 

Stt /^p 4 / n + 2 p n \ 

Introducing the dimensionless momentum x := p/m e and the parameter 6 := K(m e /Mp) n 
we have 

XF x 4 ( 1 + ( j(l!+2) a 



P = ^/ = (8) 

A white dwarf's mass depends on the star's composition. To express the white dwarf 
density then it is conventional to parameterize the density in terms of the number density 
of electrons and the atomic mass unit m u = 1.66 x 10 -24 g, so that p = pm u n e , where 
p is the electronic molecular weight. Since white dwarfs are composed mostly of carbon 
with traces of helium, the proton-electron density ratio is close to 2 since these nuclei have 
electron-nucleon ratios of 1 to 2. 

Before setting up the numerical solution, it is useful to estimate the affect of the modifi- 
cations. From equation (JSJ), one can show that the pressure in the ultra- relativistic limit is 
approximately 



3/i 3 



: 2(n + l) 
(n + 4) 

Thus, the MDR effectively add an attractive (for k < 0) or repulsive (for k > 0) force that 
scales with the Fermi momentum as p F + . 

In equilibrium the gravitational attraction must be balanced by the outward force due 
to the pressure of the degenerate gas. So, up to factors of order one, the ratio GM 2 / R A P 
should be constant. Using the pressure in the ultra-relativistic regime and p ~ M/R 3 we 
see that this ratio becomes 



2 



GM 
~R T P 



oc 



M 2/3 



M^ 3 h 
rnJ^MpR 



(10) 



neglecting factors of order one. In the non-relativistic case the ratio scales as M 1 ^ 3 R so 
increasing mass requires decreasing radius to maintain the value of the ratio. This drives 
the system into the relativistic regime given in equation (fTU|) . The leading term is constant 
in R indicating that here are equilibrium solutions for suitable mass - the mass limit - and 
radius. However, if the mass exceeds the limit then the gravitational force will exceed the 
pressure and the star will collapse. The correction due to MDR modifies this description. In 



equilibrium the ratio of equation ( TTOj) has a definite value, as before. For k > as the mass 
M is increased, equilibrium may be restored by decreasing the star's radius R, achieving a 
new equilibrium at higher density. For k < 0, as the mass is increased the correction term 
must play a role and the star must increase in relative size to remain in equilibrium. For a 
star at the Chandrasekhar mass limit the correction, approximately (M/m u ) n ^ 3 (Lp/ R) n , is 
negligible, only 10~ 14 for n = 1. (Lp is the Planck length.) We will see that this is born out 
in the numerical analysis. 

At T = all states are inside a sphere of radius pf in momentum space so the number 
density is simply 



giving a density of 



where 



-=4(t)' <»> 



87r /m e \ 3 , , 
P = -^-pm u I — J = fip x F (12) 



3 ^ "V /i 



po-=Y mu \T) - 9 - 81x10 s cm (is) 



9[), the star is in equilibrium when 



and Xp — pp/m e . 

In the context of the Chandrasekhar mass limit (see 
the gradient of the degeneracy pressure supports the Newtonian gravitational attraction. For 
a spherically symmetric mass distribution the radial pressure gradient satisfies 

dP Gp(r)m(r) 



dr r 2 
Differentiating and gathering terms gives 



(14) 



1 d fr 2 dP\ 

-— -— = -AirGp. (15) 



r 2 dr \ p dr J 

Note that we have not made the approximation that the density p(r) is uniform. 

To set up the numeric calculation it is useful to introduce the dimensionless electron 
energy z(r) = E(r)/m e and its value at the center of the star, z c := z(0). A convenient 
dimensionless radius ( is given by 

r 1 M P h 1 n 1 

(■= — with r Q := — - — ~ 7.77 x 10 9 cm. (16) 

r Q 2V3 iTLu m e pz c pz c 

Normalizing the dimensionless energy z(() to the value at the center z c gives the normalized 
energy Q(() := z(Q/z c . The radial evolution of this normalized energy is determined by the 



equilibrium condition ( TT5|) . To derive the resulting equation note that 
1 dP dx p dz 



ldP 

p dr 



pdxF dz dr 
1 



(17) 
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1 + 5 



(n+2) n +l 
2 F 



l + x 2 + 5xl +2 



x F + 6^x n F +1 



dz 
dr 



dz 



(18) 



3p a ph 3 J dr 

The third line follows from the fact we used the group velocity in the expression for the 
pressure. Using this result and the dimensionless radius ( in equation (fl5j) gives 

'x F {Q) x ° 



d 2 Q 2 dQ 







(19) 



where Xp(Q) is the dimensionless momentum in terms of Q{() as determined by the equation 

1 



Q z = -(l + x 2 F + 5x n F +2 ) 

Z c 



In the usual context of special relativity when k 



(20) 

0, xp{Q) is simply given by 



\/ Q 2 — l/^c- For large values of the central energy, it is possible to approximate the xf(Q) 
as 

/ r r 1 / . i \ n/2 ~ 

(21) 



Xf 



i Q 2 - i 



1 



- 2 K 



Q 2 



The boundary conditions for the evolution equation f fT9|) at C = are Q(0) = 1, by definition, 
and dQ/d(\^ = o = 0, to ensure vanishing pressure gradient as required by equation (JHJ). 

The non-linear differential equation f|T5|) may be solved numerically, once values for the 
parameters n, 5, and z c are chosen. To find solutions, the normalized energy Q is evolved 
via equation f[T9~|) outward from the center of the star. The outer boundary is determined by 
vanishing density, which in terms of the normalized energy means that xf{Q) = 0. When 
z c is very large, this is well approximated by Q(Cr) ~ 0- F° r an example numerical solution 
consider the case n — 1, K — —1, and z c = 10 21 . First the physical root of the cubic equation 
for xf{Q), equation (|20|) . is found by selecting the root for which the density is positive and 
small for small Q (since p cx x 2 F ~ Q 3 ). Using this root the resulting equation f[T9~|) is 
numerically evolved outward until Q{(r) — 0. The resulting numerical solution is shown in 
figure [1] as the (red) leftmost curve. The other solutions for n = 1 are found in a similar 
manner. The case for which n — 1, k — +1, and z c = 10 21 is the (green) rightmost curve. 
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The mass of the star is determined by integrating the density to the surface of the star 



r = R 



HE 



M = 4vrr 2 p(r)dr = 4npp a 3 / £ 2 x 3 F (()d( 
Jo Jo 



Anpz c p ot 



o dC, \ c(C ' 



\ ! •' /H I 

(22) 

<R 



I 1 ) ™l L 



2 dQ 



where the second line follows from equation (TTTJT) . Thus for large z c , the product of the root 
of Q(C) = an d the slope at this radius give the mass of the star. 

The mass limit M c h is determined in the ultra-relativistic limit. This is often denoted 
z c — > oo. However in the present case the onset of the MDR contributions to the pressure 
occur at central energies of z c ~ 10 16 . Since this energy is so well separated from the usual 
scale z c ~ 10 4 , we can clearly see the regime of the usual Chandrasekhar mass limit and a 
new regime where the MDR effects are manifest. 

For the k — 1 case plotted in figure [T] (n = 1, z c = 10 21 ) we find that, within the extent 
of the expected validity of the MDR, the new mass limit Mmdr becomes 

M MD r= L618 (f) 2 M & ( 23 ) 

about 10 % larger than the usual, 1.456M , result. The case for which k = — 1 and 5z c = 
—0.1 gives 

Mmdr = 1-333 2 M . (24) 

Notice that it is possible to either increase or decrease the Chandrasekhar mass limit. 

The radius of the star is found from £r, R = r (n- As in the usual case, because of the 
scaling of z c in r Q , the radii vanish as z c — > oo. In the above examples the values for Rp are 
6 x 10~ 12 cm and 8 x 10~ 12 cm, respectively, as anticipated in the qualitative argument. It 
is clear from these figures that new physics, such as the general relativistic correction (the 
Schwarzchild radius is a few kilometers), dominates well before the MDR effects become 
important. 

The density p/(pp ) = (xf(Q)/z c ) 3 for the solutions of the ±1 and k = 

(\Sz c \ = 0.1) are plotted in figure El They are manifestly non-constant. 

The solution space of equilibrium configurations is confined by the regime of validity of 
the MDR, 5z™ < 1. Within this limit the n = 1 case was solved using the exact expression 



for the physical root xf(Q) of the cubic equation. Solutions showed corrected mass limits as 
the central normalized energy rose above 10 16 . As one might expect, the corrections become 
significant as 5z c increases to 1/2, when the mass limit becomes about 1.94 solar masses for 
k > and .84 for k < 0. So as the central energy increased from ~ 10 16 the MDR effects 
continued to increase (or decrease) the mass limit to the limits of the effective description 
at onset of the Planckian regime. 

The n > 1 solution spaces displayed qualitatively similar effects on the mass limit. Using 
the approximation of equation (12TI) . the solutions were obtained for n = 2,3, 4. For example, 
the mass limit for n — 3, k — 1, z c — 3 x 10 21 becomes 1.557(yu/2) 2 M . 

In summary, the numerical solutions show that the Chandrasekhar mass limit is altered 
by about 10/%due to the modifications in the dispersion relations, increasing for positive k 
and decreasing for negative k. However, these effects occur at such high central densities 
that dynamical instabilities occur well before such energies are reached, rendering the MDR 
effects irrelevant for astrophysics. 



III. CONCLUSION 



Unlike the threshold calculations, which are rich in new features and lead to strong 
constraints on the parameters (e.g. 3j), the calculation of the Chandrasekhar mass limit 
using MDR is straightforward. The mass limit is raised or lowered according to the sign of 
the modification. As one can see in the numerical results, the qualitative argument, or from 
equation fT2Tl) for high z c , the effects become important when 5z™Q n ~ 1 or near the Planck 
scale, E/Mp ~ 1. This is well-beyond the scale at which other well-known processes occur 
so these new mass limits are not directly astrophysically relevant. 

This exact analysis corrects the results reported in pj , where it is suggested that current 
data rule out the k > case (a < in 6]). In the analysis of [6j the MDR correction to the 
radius, 5R, is 

5M P \m u 



1/3 



M ) 



2/3 



-3/2 



(25) 



where M 2/3 = (ft/37rG)(97r/8m u ) 4/3 . Hence, as M -+ M the MDR correction grows arbi- 
trarily large. This observation led Camacho to conclude that current data seems to rule out 
positive values of k. This conclusion rests from the incorrect assumption that the density is 
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FIG. 1: Results from numerical integration of the equilibrium condition for n = 1, K = ±1, and 
z c = 10 21 . Solutions for the normalized energy Q(C) are plotted for k = — 1 the lower (red) curve, 
K = - the unmodified case - the middle (blue) curve, and n = +1 the upper (green) curve. The 
inset plot shows the neighborhood of £ = —7. This plot shows the zeros of Q and the relative 
slopes, important for determining the limiting mass, for the cases k = —1, k = 0, k = +1, from left 
to right. The color coding is the same as in the larger plot. In addition, the numerical solution 
with approximation of equation (|21h is included for the k = 1 case. At this scale it is barely 
distinguishable from the solution with the exact relation for the density. 

uniform. As the present analysis shows (see figure [2]) the density is not constant throughout 
the star. While the radius does increase with central energies for positive k, the actual radii 
at these energies are so small that other processes such as neutronization and gravitational 
collapse occur long before the star evolves to that state. 
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FIG. 2: The densities p/(pp a z^) are shown for \Sz c \ = 0.1. On the left of the plot the solutions 
are, from top to bottom k = —1 (red), n = (blue), and k = +1 (green). 



Current observations seem to indicate white dwarfs with smaller radii than expected, at 



least for iron cores 



10 



ll|. The conclusion of [6] suggests that solutions with negative k 



(or a > 0), when the white dwarf radius in the MDR case is reduced, lead to a possible 
explanation for the smaller radii. However as seen in the result reported after equation (1241) . 
the model with the realistic density shows significant corrections only at the Planck scale. 

The argument rests on fermion statistics and the modified dispersion relation. Hence, in 
so far as these relations hold, it would seem that the calculations of the mass limits generalize 
to frameworks outside the EFT approach. However, as pointed out in 16[, this neglects the 
possible changes in the integration measure of the pressure, equation [7J For instance, in 
deformed special relativity (DSR) [3, [l^ the modified dispersion relations are regarded as 
the invariant of a relativity group with two invariant scales - not only the speed of massless 
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modes, c, but also an invariant length or energy 19[ . As interpreted by Hossenfelder 15] the 
DSR dispersion relations arise from an effective description of gravitational effects at high 
energy- density. The calculation of the neutron star mass in DSR was carried out in [bsl ]. 

So barring modification of statistics or of the measure, the present calculations are valid 
mutatis mutandis in a more general context. The new mass limits would appear not to 
have astrophysical relevance. Neutronization and the effects of general relativity would still 
occur at much lower momenta. Recently there has been work on the impact of MDRs on 
the structure of astrophysical objets [l7]. 

The MDR effects on the Chandrasekhar mass limit may also be interpreted in the context 
of generalized uncertainty principle. The non-linear relation between particle momentum 
in a high energy- density region and the asymptotic region, induces a modification in the 
Heisenberg uncertainty principle 15] 

h 



AxAp > 



dp 
dk 



(26) 



Expanding the modification as dp/dk = 1 + n(p/Mp) n and following the method of estimat- 
ing the ground state energy of the hydrogen atom from the uncertainty principle, one may 
show that the pressure is modified as in equation (J7J) (up to a factor of 1/2). This leads to 
the same effects as derived in the EFT context. 
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